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CHAPTER  I 


PREFACE 

This  work  attempts  a  broad  analysis  of  the  form  of 
pupillary  response  to  light  and  the  mechanisms  which  combine  to 
determine  this  response.  The  literature  is  voluminous,  though 
quantitative  data  had  to  await  the  development  of  suitable  pupil- 
lometers  in  more  recent  times.  Many  studies  have  been  extensive  in 
both  the  experimental  and  theoretical  realms.  It  is  surprising,  in 
fact,  the  amount  of  interest  in  this  system  which  has  captured  the 
minds  of  scientists  throughout  the  ages.  Modern  times  has  seen  some 
very  fine  and  extensive  experimental  investigations  as  well  as  some 
well  worked  out  theoretical  formulations.  This  work  will  combine 
these  two  approaches. 

The  pupil  system  possesses  certain  properties  which  make  it 
attractive  for  a  study  of  neural  reflex  and  control. 

1)  The  end  organs  are  superficial  and  readily  available 
without  the  need  for  dissection. 

2)  Though  the  retina  is  acknowledged  to  be  a  highly 
elaborate  structure,  the  small  number  of  central 
synapses  make  this  relatively  uncomplicated  as  a  CNS 
system. 

3)  Although  secondary  factors  such  as  fatigue  and  alert¬ 
ness  of  the  subject  influence  pupillary  activity,  the 
fact  that  learning  apparently  does  not  occur  simpli¬ 
fies  analysis  of  data  and  experimental  design  using 
human  subjects. 
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It)  As  opposed  to  many  other  reflexes  this  has  a 
readily  defined  teleological  foundation. 

Thus,  the  system  is  simple,  functionally  complete  and 
readily  manipulatable,  and  as  has  been  mentioned,  investigators 
have  not  neglected  these  virtues. 

We  will  precede,  then,  with  a  review  of  some  of  the  perti¬ 
nent  anatomy  and  physiology  and  present  some  new  experimental  work. 
Using  concepts  from  systems  analysis,  we  will  then  investigate 
theories  of  system  operation  under  certain  circumstances  and  draw 
some  broad  conclusions. 
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CHAPTER  II 


ANATOMY  OF  THE  SYSTEM. 

Parasympathetic  system: 

The  pupillary  light  reflex  encloses  an  arc  from  retina 
through  midbrain  to  pupil.  Retinal-  afferents  probably  perform  the 
dual  role  of  conveying  visual  information  and  carrying  light  in¬ 
tensity  information  for  the  pupil  response  (Schweitzer,  1956). 

Optic  nerve  fibers  send  off  collaterals  which  project  to  the  region 
of  the  rostral  portion  of  the  third  nerve  nucleus  (Ranson  and 
Magoun,  1933) •  Many  investigators  feel  that  there  is  a  single 
synapse  in  the  pretectal  region,  though  Polyak  (1957)  presents 
evidence  that  the  pregeniculate  nucleus,  anatomically  closely  as¬ 
sociated  with  the  main  geniculate  nucleus,  acts  as  a  way  station  for 
the  pupillary  tract  and  sends  axons  to  the  occulomotor  nucleus  by 
way  of  the  pretectal  region  without  synapse.  In  any  event,  there  is 
agreement  that  there  is  decussation  of  the  pathways  from  the  two 
geniculate  bodies  prior  to  synapse  on  motor  ganglion  cells.  The 
parasympathetic  final  common  pathway  therefore  receives  contra  and 
ipsilateral  innervation  at  this  point,  and  the  consensual  reflex  is 
established.  Pupil  activity  on  the  two  sides  is  highly  correlated. 
Cells  in  the  rostral  part  of  the  third  nerve  nucleus,  generally 
thought  to  be  located  primarily  in  the  well  delineated  Edinger- 
Westphal  nucleus,  send  preganglionic  parasympathetic  motor  fibers  to 
the  ciliary  ganglion  by  way  of  the  third  nerve.  After  a  synapse  in 
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that  structure,  the  reflex  pathway  is  completed  through  the  short 
ciliary  nerves  which  pierce  the  sclera  at  the  posterior  pole  of  the 
globe  and  travel  to  the  sphincter  ptqpilli. 

Sympathetic  system; 

Without  doubt  the  sympathetic  nervous  system,  plays  a  part 
in  determining  level  of  operation  and  dynamic  behavior  of  the  pupil 
under  certain  conditions,  as  has  been  substantially  shown  by  Lowen- 
feld  (1958).  She  has  brought  together  evidence  indicating  the 
presence  of  a  psychosensory  center  in  the  hypothalamus,  and  she  has 
placed  on  a  firm  foundation  the  presence  of  a  discrete  radial 
dilatory  muscle.  The  outflow  of  the  sympathetic  system  proceeds 
mainly  through  the  first  two  thoracic  segments  up  the  chain  to 
synapse  in  the  superior  cervical  ganglion.  Postganglionic  fibers 
travel  into  the  cranium  on  the  carotid,  leaving  the  artery  before 
its  termination  to  seek  a  path  crossing  the  floor  of  the  middle  ear 
and  thence  to  join  the  third  nerve  distal  to  the  gasserian  ganglion. 
Long  ciliary  nerves  carry  the  signal  to  the  globe  partly  bypassing 
and  partly  through  the  ciliary  ganglion  without  synapse. 
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CHAPTER  III 


REVIEW  OF  THE  PHYSIOLOGY 

Pupillary  noise: 

Under  average  lighting  conditions  and  fixed  accommodation, 
one  notices  a  persistent  fluctuation  in  the  size  of  the  pupil  not 
obviously  correlated  with  other  physiological  activity.  Stark, 
Campbell  and  Atwood  (1958)  have  pointed  out  that  there  is  high  cor¬ 
relation  between  the  fluctuation  of  the  two  pupils.  Noise  power 
increases  with  retinal  illumination.  The  spectra  of  the  noise  signal 
and  pupil  response  have  different  break  frequencies  with  parallel 
attenuation  and  therefore  require  partially  separate  channels.  Con¬ 
clusions  as  to  where  the  noise  could  be  injected  into  the  system 
and  the  nature  of  the  source  were  reviewed.  Lower  stein  (1950)  has 
shown  much  decreased  noise  activity  after  lesions  in  the  hypothalamic 
center  of  pupil  sympathetic  activity  in  monkeys,  implicating  the 
sympathetic  as  the  more  probable  pathway. 

Pupil  threshold: 

Work  on  pupil  threshold  to  increases  of  light  intensity 
has  been  reviewed  by  Schweitzer  (1956),  who  also  contributed  experi¬ 
mental  evidence  carried  out  with  an  infra  red  optical  system,  deter¬ 
mining  threshold  activity  by  subjective  methods.  For  small  simulating 
fields  pupil  threshold  was  found  to  be  much  below  visual  threshold, 
but  as  the  field  was  enlarged  sensitivity  of  the  pupil  system  in¬ 
creased  and  with  large  fields  approximated  visual  thresholds. 
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Therefore,  unlike  the  visual  system,  spatial  summation  is  a  funda¬ 
mental  property  of  the  pupil  system.  As  to  temporal  summation, 
though,  the  two  systems  act  similarly.  Within  a  specific,  small 
interval,  measured  by  Bouman  (1950)  to  be  0.02  and  0.035  secs, 
cases  threshold  depends  on  total  energy  independent  of  time  distri¬ 
bution,  (the  Bunsen-Roscoe  Law).  As  the  time  distribution  of  energy 
increases  beyond  that  critical  value  threshold  energy  values  in¬ 
crease.  These  two  insights  into  threshold  dependence  on  time  and 
area  parameters,  reveals  the  fact  that  constant  energy  stimuli  yield 
constant  responses,  given  certain  bounding  limits  of  stimulus  duration 
and  space  variables.  Threshold  experiments  have  also  revealed  that 
the  pupillary  reflex  is  not  dependent  on  one  class  of  receptor 
(Schweitzer,  ibid;  Lowenstein  and  Lowenfeld,  1959)?  though  Lowen stein 
notes  that  the  red,  as  opposed  to  the  white  light  gives  the  larger, 
more  sustained  response. 

Response  to  static  input: 

When  light  is  suddenly  applied  to  the  dark  adapted  eye  and 
then  maintained  constant  for  a  period  of  time,  there  results  an 
initial  sharp  constriction  which  may  be  maintained  in  a  plateau  for 
several  seconds,  subsequently  followed  by  a  slow  dilation  over  a 
period  of  minutes  to  some  relatively  fixed  area,  less  than  the  dark 
adapted  area.  There  is  great  variability  in  the  data  among  various 
investigators  of  the  static  response  (reviewed  by  Schweitzer),  but  in 
general  plots  of  pupil  area  vs.  wide  field  illumination  result  in 
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sigmoid  shaped  curves.  (Reeves,  1918).  The  teleological  effect  of 
such  a  pattern  of  response  is  to  allow  effective  functioning  over  a 
very  broad  operating  range. 

Dynamic  response  -  steps  and  impulses? 

Talbot  (1939)  used  a  large  field  for  stimulation  and  re¬ 
corded  pupil  diameter  photographically.  He  described  a  response  to 
a  step  input  of  light  which  began  after  a  latency  of  0.25  secs., 
occasionally  preceeded  by  a  small  dilation,  showed  a  '‘jog1*,  or  rapid 
change  of  rate  of  contraction  after  0.6  secs,  with  decreased  velocity, 
and  a  second  jog  at  0.9  secs.  The  transient  was  completed  in  2.5  secs. 
For  small  responses,  less  than  a  1.5  mm.  change,  the  peak  was  followed 
by  an  immediate,  relatively  rapid  dilation  or  rebound.  Latency  of 
the  dilatory  phase  following  step  decrease  of  light  was  O.h  secs, 
at  an  initial  speed  of  0.25  mm/sec.,  half  the  rate  of  initial  con¬ 
striction.  The  response  to  an  impulse  of  light  showed  the  latency, 
speed,  and  the  characteristic  discontinuities  seen  in  the  step  response. 
Of  particular  note  was  the  observation  that  the  minimum  duration  of 
the  transient  was  1.1  secs.  For  short  flashes  (up  to  1  sec.  duration), 
he  found  that  the  amplitude  of  the  response  varied  with  the  stimulus 
energy,  but  that  there  was  a  falling  off  of  response  at  higher  in¬ 
tensities.  In  an  elaborate  theoretical  discussion  he  postulated  an 
interplay  of  the  two  muscles,  radial  and  constrictor,  which  had 
mechanical  properties  of  tension  and  elasticity.  A  proper  coordina¬ 
tion  of  excitation  times  would  then  give  rise  to  the  characteristic 
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jogs,  which  played  an  important  part  in  his  scheme.  He  also  as¬ 
sumed  that  the  frequency  of  nerve  impulses  was  proportional  to  log 
intensity  with  rapid  formation  of  chemical  muscle  excitant  and 
therein  derived  an  exponential  formula  for  contraction  with  a  time 
constant  dependent  on  step  intensity. 

Van  der  Tweel  (1959)  was  also  interested  in  developing  a 
model  of  the  pupil  system  based  on  experimental  data  of  the  overall 
response  to  steps  and  impulses  and  assumptions  about  the  behavior 
of  the  system  elements  inferred  from  other  general  physiological 
work.  His  experimental  method  used  an  optical,  system  in  which  the 
pupil  was  illuminated  by  deep  blue  light  and  the  image  of  diameter 
was  projected  onto  moving  blue  sensitive  film,  thus  giving  a  con¬ 
tinuous  record  of  pupil  diameter.  Unfortunately  the  resolution  of 
the  technique  was  poor,  but  his  data  is  somewhat  similar  to  Talbot* s 
and  ours  on  points  of  comparison.  His  data  seemed  to  indicate  that 
the  peak  of  responses  to  impulse  stimuli  with  durations  up  to  0.5 
sec.  had  a  constant  time  duration  to  peak,  and  he  thought  that  within 
this  time  range  the  stimulus  appeared  to  be  a  trigger  mechanism, 
the  height  of  the  response  being  related  to  the  initial  flux  by  a 
log  function.  A  second  flash  applied  0.3  -  0.5  secs,  after  the 
first  showed  summation  logarythmically.  Also,  he  reasoned  from  the 
fact  that  since  time  constants  of  contraction  to  both  large  and 
small  stimuli  was  constant,  the  motor  element  could  not  be  amplitude 
limited  in  the  range  he  investigated.  Unlike  Talbot  and  Lowenstein 
(see  below),  in  considering  dilation,  van  der  Tweel  assumed  that 
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dilation  was  passive,  based  mainly  on  the  fact  that  his  data 
demonstrated  that  the  speed  of  dilation  depended  on  the  peak 
diameter  reached  following  an  impulse.  He  noted  that  dilation 
approximately  followed  an  exponential  curve, 

Lowenstein  and  co-workers  have  done  much  work  on  phases 
of  the  pupil  over  a  period  of  many  years  using  cats,  rodents,  birds 
and  primates;  and  neurophysiological  techniques  of  injury  ablation 
and  pharmacological  devices.  His  instrumentation  has  included  infra¬ 
red  photography  and,  recently,  an  infra-red  scanner  which  selects 
the  pupil  diameter  from  a  sampling  of  pupil  segments  made  ten  times 
a  second.  He  has  studied  the  form  of  the  response  with  a  view 
to  the  fine  points  of  the  pattern  and  has  given  much  weight  to  the 
first  derivative  of  the  response  curve  as  an  indicator  of  fundamental 
physiological  processes.  The  light  reflex  in  his  conception  is 
finely  balanced  between  sympathetic  and  parasympathetic  control 
(Lowenstein  and  Lowenfeld,  19!?0);  interference  in  this  balance  pro¬ 
duces  characteristic  changes  in  reflex  pattern.  In  the  normal  pupil 
contraction  is  in  a  single  wave  of  acceleration  and  deceleration  to 
a  peak  occurring  several  tenths  of  a  second  after  the  cessation  of 
the  one  second  stimulus  impulse.  Redilation  is  achieved  in  two 
phases.  The  first  is  fast  and  is  primarily  parasympathetic  inhibi¬ 
tion.  The  second  is  slower  and  requires  an  intact  sympathetic  con¬ 
nection  for  its  completion.  Peripheral  sympathectomy  has  an  effect  on 
the  contraction  phase,  viz.  1)  latency  of  contraction  is  shorter; 

2)  amount  of  contraction  is  greater;  3)  speed  of  contraction  is  greater; 
and  U)  time  to  peak  is  shorter.  An  interesting  change  in  the  response 
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occurred  after  hypothalamic  lesions  in  monkeys.  The  rate  of  con¬ 
traction  was  increased,  but  constriction  halted  abruptly  before  the 
end  of  the  one  second  stimulus.  The  pupil  maintained  contraction  on 
a  plateau  until  0.1  or  0.2  tenths  seconds  after  the  end  of  the 
stimulus,  and  then  redilation  took  place  rapidly.  Lowenstein  and 
Lowenfeld  (1952)  have  correlated  changes  in  the  response  with  fatigue 
with  responses  after  CNS  lesions  and  have  thereby  formulated  an 
hypothesis  of  the  temporal  development  of  fatigue:  central  structures 
fatigue  before  peripheral  elements  and  the  sympathetic  system  before 
par a sympathetic.  Psychosensory  events  such  as  a  loud  noise  will 
restore  the  failing  reflex. 

Dynamic  response  -  sinusoidal  input: 

Stark  and  Sherman  (1957)  initiated  a  study  which  has  become 
the  most  thorough  investigation  of  the  sinusoidal  response  of  the 
pupil.  They  reasoned  that  since  the  pupil  controlled  retinal  flux 
which  was  in  turn  a  factor  in  pupil  response,  the  dynamic  operation 
of  the  system  should  most  properly  be  studied  with  the  aid  of  servo- 
analysis,  the  science  of  feedback  control  systems. 

A  review  of  the  properties  of  a  servomechanism  is  in  order 
here.  An  arch-typical  servosystem  is  the  thermostatically  controlled 
heating  system  of  a  house  in  which  a  device  samples  room  temperature 
and  operates  a  second  device  which  adjusts  room  heat  so  that  the 
activation  of  the  sampling  device  is  minimal.  This  may  be  represented 
schematically  as  shown  in  Fig.  1.  The  boxes  represent  elements  in 
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Fig-  1.  Simple  servomechanism.  Output  quantity  compared 
with  desired  reference  value  and  difference  used 
to  drive  output  device.  Arrows  designate  uni¬ 
directional  signal  flow. 
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the  system  and  the  connecting  lines  with  arrows  denotes  the  path 
of  the  signal,  or  information.  The  left  hand  junction  of  the  three 
line  segments  at  the  circle  enclosing  the  minus  sign  represents  a 
subtraction  operation  taking  place  between  the  reference  value  and 
the  sampled  value  of  the  controlled  quantity.  Therefore,  in  our 
representative  control  system,  temperature  is  the  output  quantity. 

The  thermostat  samples  the  temper  at  tare  of  the  room  and  compares 
the  resulting  value  with  a  predetermined  desired  value,  the  reference. 
If  an  error  results,  the  output  device  is  activated  to  adjust  the 
temperature  to  a  new  value  by  putting  in  or  removing  heat.  In  recent 
years  the  theory  of  such  devices  has  been  worked  out  to  a  high 
degree  for  the  linear  case  (for  example.  Bower  and  Schultheiss, 

1958s  and  Chestnut  and  Mayer,  1951).  Important  properties  are  speed 
of  response,  minimization  of  error,  and  stability. 

In  applying  servo theory  to  the  pupil  system.  Stark  and 
Sherman  were  confronted  with  the  problem  of  describing  the  control 
process  which  was  taking  place  at  the  plane  of  the  pupil.  If  a  broad 
source  projects  a  given  intensity  of  light  onto  the  whole  eye,  the 
flux  which  falls  onto  the  pupil  will  be 

F  s  Axl 

in  obvious  symbolism.  Any  change  in  F  may  be  described  by  a  change 
in  A  or  a  change  in  I.  Therefore,  in  difference  notation: 


AF  s  IxAA  =  AxAI 

If 

AF^s  AxAI 

and 

A  F0»  IxAA 

AFrAFi+AF0  o 


then 
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A  step  change  of  light  due  to  a  change  in  I  would  travel  around  the 
control  loop  and  cause  a  change  in  A.  The  above  notation  allows  us  to 
calculate  the  gain  as 

AF0  IxM 

°  ~  AxZl” 

This  formulation  has  the  difficulty  that  one  is  forced  to  choose  what 
will  be  the  intensity  I  to  be  used  in  calculating  the  numerator,  i.e., 
the  I  before  the  change,  or  after  the  change.  Those  authors  worked 
mainly  with  sinusoidal  inputs  which  varied  about  a  mean  T,  If  the  peak 
to  peak  change  (  I)  could  be  kept  small,  the  error  introduced  by 
this  complication  would  be  negligible. 

The  initial  work,  then,  described  the  properties  of  the 
open  loop  system  response.  The  data  indicated  an  amplitude  response 
break  frequency  at  1,5  cps  and  a  noruninimum  phase  shift  which  sug¬ 
gested  asymptotic  limiting  at  -270°,  Justifying  a  linear  system  on 
their  small  signal  approximation,  they  deduced  from  servo theory  a 
model  system  which  had  three  lags  and  a  gain  of  0.16,  Their  phase 
data  gave  a  transport  delay  of  0,18  secs.,  and  after  postulating 
equal  time  constants  for  all  lags,  arrived  at  the  open  loop  transfer 
function 

G(S)  -  0,l6e~0*^0 
(1+0.  IS)3" 

This  is  merely  the  Laplace  transform  of  the  system*s  response  to 
sinusoidal  inputs.  The  transform  contains  all  the  necessary  para¬ 
meters  for  completely  describing  the  behavior  of  a  linear  system,  and 
therefore  by  proper  manipulation  will  give  the  differential  equations 
of  response  in  time  (Stark,  1959). 
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Stark  and  Baker  (195>8) ,  using  arguments  from  servo theory 
and  employing  drugs  to  shift  system  parameters  explained  the  origin 
of  the  previously  described  "induced  hippus",  the  l.ii  cps  occillation 
of  the  pupil  seen  under  special  lighting  conditions. 

Tyler,  working  in  this  laboratory,  critically  reviewed  the 
problem  of  representing  the  pupil  system  by  a  linear  lumped  parameter 
network.  Working  with  his  own  data  on  sinusoidal  and  step  on  and  off 
response,  he  concluded  that  there  was  no  real  asympotote  in  the  phase - 
frequency  curve  and  had  difficulty  correlating  frequency  response  with 
step  response  data.  Such  conceptual  complications  led  him  to  con¬ 
struct  a  model  on  feed-forward  lines.  Though  his  objections  were 
aptly  presented,  his  discussion  was  along  highly  theoretical  lines 
and  without  real  physiological  basis. 
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CE4PTER  IV 


PRESENT  APPROACH 


From  the  foregoing  discussion  we  have  seen  that  there  are 
several  factors  which  determine  the  area  of  the  pupil: 

1)  Intensity  of  light  falling  on  the  retina,, 

2)  area  of  retina  illuminated* 

3)  past  and  present  time  course  of  the  light  intensity 
at  the  retina* 

10  state  of  adaptation  of  the  retina* 

3)  distribution  of  the  light  flux  between  the  two  retinas* 
6)  internal  random  influence  on  pupil  area,  or  noise. 

As  we  have  seen  from  Schweitzer 8 s  work*  1)  and  2)  may  be 
combined  as 


F  =  A*  x  I* 


where  I*  ®  intensity  of  light  at  the  plane  of  the  retina  and  A*  s 
the  area  of  retina  illuminated.  We  will  not  consider  the  effect  of 
binocular  imbalance*  point  5),  which  is  too  complicated  for  this 
preliminary  study  and  may  be  controlled  by  standardizing  experimental 
conditions.  In  point  6)  we  recognize  the  fact  that  other  influences 
than  light  level  have  a  part  in  the  determination  of  pupil  area. 
Accommodation  is  one  of  these  forces  and  may  be  controlled  experi¬ 
mentally;  fatigue  and  psychosensory  phenomenon  can  be  controlled  less 
well;  and  we  are  finally  left  with  ’‘noise",  the  random  uncorrelated 
variations  in  pupil  area  and  a  subject  for  study  in  itself.  Finally, 
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since  the  state  of  retinal  adaptation  is  a  function  of  the  time 
course  of  incident  flux,  we  see  that  pupil  area  at  any  moment  in  a 
controlled  situation  is  ideally  determined  only  by  the  past  and 
present  course  of  incident  flux,  or  symbolically : 

A“  f(F) 


The  fundamental  problem: 

The  simplicity  of  the  foregoing  conclusion  draws  us  to 
the  brink  of  the  total  problem  of  the  physiology  of  the  pupillary 
reflex  to  light,  almost  by  a  process  of  reductic  ad  absurdum, 
namely:  how  can  we  determine  the  response  to  light  at  any  time  given 
any  previous  light  flux  history?  The  answer  to  such  a  question  would 
presumably  lead  through  all  the  ramifications  of  submicro sc opic 
structure  and  function.  It  is  to  the  delineation  of  the  nature  of 
the  problems  raised  in  a  search  for  a  general  theory  that  we  are 
devoted  here. 


An  apology: 

The  central  reason  for  the  rise  of  mathematics  in  any 

science  is  no  better  stated  than  by  Philip  E.  B.  Jourdain  (195>6) : 

•’Our  ideal  in  natural  science  is  to  build  up  a 
working  model  of  the  universe  of  the  sort  of 
ideas  that  all  people  carry  about  with  them 
everywhere  *in  their  heads, 1  as  we  say,  and  to 
which  idea  we  appeal  when  we  try  to  teach  mathe¬ 
matics.  These  ideas  are  those  of  number,  order, 
the  numerical  measures  of  times  and  distances,  and 
so  on.n 
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In  short,  even  if  a  mathematical  construct  adds  nothing  new  to  our 
concepts,  such  a  representation  is  valuable  for  its  clarity,  economy, 
generality,  and  systemization,  for  it  is  in  such  an.  atmosphere  that 
correlations,  the  reasons  for  science,  are  made. 

Biology  has,  to  a  great  extent,  resisted  the  inroads  mathe¬ 
matics  has  made  into  the  physical  sciences.  Much  of  the  tardiness 
has  undoubtedly  been  due  to  the  complexity  of  the  material  and  its 
refractoriness  to  experimental  control,  but  as  Margeneau  says  (1950): 

’’Talk  about  the  complexity  of  organisms,  about  the 
sum  that  is  more  than  its  parts,  usually  winding  up 
in  the  reminder  that  biology  is  a  field  for  which 
casual  analysis  is  inadequate,  proves  nothing. 

Aristotelian  physicists  have  said  the  same  thing 
about  inorganic  nature.  Indeed  every  field  of  human 
interest  in  which  significant  variables  (variables 
of  state)  have  not  been  discovered  is  for  that  very 
reason  complex  and  intractable,  and  to  say  that  it 
defies  casual  analysis  is  merely  to  admit  that  it 
has  thus  far  not  been  penetrated  by  the  usual  methods 
of  science.  ” 

In  studying  the  pupillary  reflex  to  light  we  are  dealing 
with  a  system  which  is  acted  upon  by  a  forcing  variable,  light,  and 
reacts  in  a  characteristic  response  mediated  by  a  system  of  trans¬ 
ferences  and  transformations  of  information.  Solutions  to  the 
problems  presented  by  such  a  study  parallel  in  many  ways  methods  in 
systems  analysis  and  the  science  of  response  of  physical  systems. 

Indeed,  Norbert  Wiener,  who  is  frequently  accused  of  trying  to 
fabricate  brains  out  of  machines,  in  reality  has  appealed  to  the 
economy  of  human  effort,  when  he  explains  the  touchstone  of  his  philoso¬ 
phy  in  the  Introduction  to  Cybernetics: 
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’There  are  fields  of  scientific  work,  as  we 
shall  see  in  the  body  of  this  book,  which  have 
been  explored  from  the  different  sides  of  pure 
mathematics,  statistics,  electrical-  engineering, 
and  neurophysiology;  in  which  every  single  notion 
receives  a  separate  name  from  each  group;  and  in 
which  important  work  has  been  triplicated  or  quad¬ 
ruplicated;  while  still  other  important  work  is 
delayed  by  the  unavailability  in  one  field  of  re¬ 
sults  that  may  have  already  become  classical  in  the 
next  field.” 

In  the  course  of  this  study  we  will  therefore  not  hesitate 
to  draw  from  these  fields  for  methods  of  symbolic  representation,  ex¬ 
perimental  approach,  and  formal  analysis. 


Systems  analysis  -  an  introduction. 

The  case  in  point  is  a  system  made  up  of  real  components 
strung  together  in  some  manner,  which  reacts  in  a  measurable  way  (the 
output)  to  the  stress  of  a  forcing  function  (the  input).  On  p. 
we  used  the  simple  notation  of  connecting  the  diagrams tic  representa¬ 
tion  of  two  measurable  points  of  our  real  system,  viz,  input  and 
output,  by  lines  and  boxes.  The  boxes  represented  the  real  lumped 
elements  of  the  system  which  generate,  store,  dissipate  or  transform 
energy,  or  information.  The  lines  connecting  the  boxes  represent  uni¬ 
directional  transfer  of  information  unchanged  from  one  element  to 
another.  A  simple  system  might  therefore  be  diagrammed  as  in  Fig.  2. 
where  the  input  function  is  e-j_,  which  is  transformed  into  an  output  e0 
through  two  operators  and  G2.  By  knowing  that  there  exists  two 
distinct  system  elements  we  infer,  or  perhaps  have  measured,  an  inter¬ 
mediate  quantity  e^.  Obviously,  if  we  knew  only  the  input  and  output 
of  the  system  and  nothing  of  its  mechanism  of  operation,  we  would  have 
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Fig.  2.  Two  element  series  network.  Signal  flow 
representation  of  transformation  of  e^ 
into  eQby  operators  and  G^. 
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simply  represented  the  transformation  e^->>e0  by  a  single  block,  G. 

One  can  see,  then,  that  this  notation,  the  block  diagram,  can  be 
expanded  as  our  insight  into  mechanisms  of  the  system  increases. 

Concentrating  on  the  first  block,  we  will  call  the  operator 
Gp  a  transfer  function.  Using  the  notation  of  modern  algebra,  the 
transformation  ej_— *.e0  becomes 

el  =  °1  ei 

In  general,  these  functions  are  not  commutable,  i.e. 

ei  *  ei  Gi 

in  the  general  case.  This  becomes  evident  on  the  consideration  of 
nonlinear  operators  to  be  discussed  shortly.  The  transformation  may 
also  be  visualized  in  the  usual  manner  using  the  graph  of  the  function. 


Finally,  since  eQ  “Gpe-^  , 

The  complete  system  is  represented  algebraically  as 

e  =G0Gn  e-,  , 
o  2  1  1  * 

again  noting  that  in  general 

G2Glei  ^  GlG2ei  • 

The  properties  of  the  operators  fall  into  general  classes, 
two  of  the  most  important  being  linear  and  nonlinear  operators.  Since 
there  are  well  developed  methods  in  linear  analysis  and  nonlinear 
analysis  is  still  troublesome,  difficult  or  impossible,  attempts  are 
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frequently  made  to  approximate  the  properties  of  systems  by  linear 
functions.  The  property  of  linearity  is  based  on  the  principle  of 
si^perposition  to  which  a  system  conforms  if  the  sum  of  two  or  more 
input  functions  yields  an  output  function  equivalent  to  the  sum  of 
outputs  to  the  individual  inputs,  or 
If  71  -  f(x1) 

y2  =  f(x2^ 

•  •  • 

*  0  0 

Then  y1+y2  +  .#.  rf(x-L+x2  +  ...  ) 

With  such  a  system  if  the  input  is  doubled,  the  output  is  doubled. 

In  passing  it  should  be  noted  that  if  an  operator  varies  with  time,  the 
system  may  still  be  linear  (see  Ku  for  a  discussion). 

While  for  some  purposes,  particularly  when  there  is  a  need 
for  mathematical  analysis  in  verification  or  synthesis,  one  can  throw 
away  much  important  information  when  the  linear  approximation  is  at¬ 
tempted,  in  fact  a  realization  of  the  properties  conferred  by  nonlinear 
elements  can  contribute  much  to  an  understanding  or  analytical  dissection. 

We  are  now  in  a  position  to  note  the  validity  of  the  state¬ 
ment  that  the  operators  in  the  algebraic  notation  are  not  commutable  in 
the  general  case.  If,  for  example  Cq_  was  linear  and  G2  generated  the 
log  of  its  argument,  then 

e^  =G2Gi6j_  r  log  (ae^+c) 
and  ej_  sGiG2e-}_  =  a  log  ej+c 

These  are  two  wholely  different  functional  relationships. 

Selin  (1958)  suggested  the  use  of  just  such  a  differential 

point  in  the  study  of  a  nonlinear  property  in  the  pupil  system.  He 
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reasoned  that,  for  the  purposes  of  his  analysis  the  pupil  system 
could  be  broken  down  into  two  elements:  1)  A  frequency  dependent  element, 
in  which  the  output  for  any  given  amplitude  input  was  dependent  on 
frequency,  and  2)  an  amplitude  dependent  element  which  saturated  at 
higher  values  of  input  but  was  independent  of  frequency.  These 
functions  are  represented  in  figures  3  and  U*  If  the  G2  element  pre¬ 
cedes  the  G 1  element,  the  break  point  of  the  system  transfer  function 
would  be  independent  of  the  input,  e^.  On  the  other  hand  if  the  fre¬ 
quency  dependent  component  comes  first  and  high,  saturating  levels  of 
input  are  being  used  the  output  at  low  frequencies  would  saturate  the 
nonlinear  element  while  at  high  frequencies  the  input  to  G2  would  be 
less  than  the  critical  value  ec.  The  comparison  of  the  frequency  re¬ 
sponses  is  shown  in  Fig.  3* 

A  purely  analytic  method  has  been  used  to  * dissect*  the 

system. 

The  approach: 

The  remainder  of  this  work  will  be  an  exploration  of  the 
ideas  of  systems  analysis  in  their  application  to  the  pupil  system. 

The  experimental  problems  to  be  confronted  are  the  deter¬ 
mination  of  the  functional,  relationships  governing  the  response  of 
the  pupil  under  different  operating  conditions.  Explicitly  the  ques¬ 
tions  are: 

1)  How  does  the  response  to  an  impulse  of  light  depend  on 
(a)  the  duration  of  the  input  function  and  (b)  the 
height  of  the  impulse? 
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2)  How  does  the  response  to  sinusoidal  varying  light 
depend  on: 

(a)  the  peak  to  peak  variation,  and  (b)  the  mean  value 
of  the  stimulus? 

With  this  data  at  hand  it  will  be  in  order  to  draw  some  in¬ 
ferences  about  the  operation  of  the  system. 
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Fig.  3.  Diagramatic  frequency  response  function.  Solid 
curve  represents  response  of  frequency  dependent  operator. 
Dashed  curve  represents  response  of  saturating  system 
in  which  saturating  operator  follows  frequency  dependent 
operator. 


Fig.  4.  Input-output  relationship  of  saturating  element. 
Gain  independent  of  input  for  all  values  less  than  q  ’(c)* 


CHAPTER  V 


METHODS 

Instrumentation  -  the  pupillometer : 

The  methods  used  for  stimulating  and  recording  pupil  area 
have  been  reviewed  elsewhere  (Stark,  19^9).  They  will  be  discussed 
here  for  completeness. 

The  pupillometer  is  composed  of  three  subunits,  each  with 
a  partly  independent  light  path,  but  all  ultimately  directed  through  ^ 
an  objective  lens  before  the  eye. 

1)  The  infrared  channel: 

A  beam  of  light  from  a  high  intensity  tungsten  source  passed 
through  a  infrared  filter  and  was  focused  on  the  iris  in  a  circle  always 
larger  than  the  pupil.  The  amount  of  the  infrared  light  reflected  from 
the  iridopupillary  plane  was,  to  a  good  approximation,  inversely  pro¬ 
portional  to  the  pupil  area,  and  the  intensity  of  the  reflected  com¬ 
ponent  was  detected  by  a  properly  oriented  infrared  photocell.  The 
resulting  output  was  recorded  on  a  Sandborn  pen  recorder  simultaneously 
with  a  record  of  the  light  input.  Increases  in  light  level  and  de¬ 
creases  in  pupil  area  were  recorded  as  positive  diflections,  whence 
the  terminology  ’’rising”  and  ’’falling”  phases,  etc.  The  infrared  , 
filter  in  the  impulse  experiments  was  composed  of  two  Kodak  87c  filters 
tandem.  A  faint  dark  red  could  be  seen  by  the  dark  adapted  subject 
when  the  stimulating  light  was  covered,  but  this  was  all  but  masked 
by  even  the  absolute  minimum  intensity  of  the  stimulating  light.  In 
the  sinusoidal  experiments  an  army  snooperscope  filter  was  used  and 
was  absolutely  opaque  to  visual  wave  lengths. 
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2)  Stimulating  channel: 

A  S$yl  vania  glow  modulator  lamp  was  used  as  the  stimulating 
source.  It  had  frequency  response  characteristics  well  beyond  the 
pupil  range,  and  because  its  output  could  be  completely  controlled 
electronically,  it  was  ideally  suited  for  this  application.  Its  one 
disadvantage  was  the  range  of  flux  variation  which  was  limited  by  a 
minimum  low  level  as  well  as  the  maximum.  Kodak  neutral  density 
filters  could  be  inserted  into  this  path  to  reduce  transmittance. 

The  flux  from  this  source  was  focused  at  the  plane  of  the  pupil  as  a 
small  circle  always  smaller  than  the  smallest  area  of  the  pupil.  In 
this  way  change  in  pupil  size  would  have  no  effect  on  the  amount  of 
light  reaching  the  retina,  and  this  manner  of  operation  is  called 
"open  loop1’. 

3)  Calibration  channel: 

Calibrations  used  to  correlate  absolute  pupil  area  with  in¬ 
frared  photocell  output  were  obtained  photographically  along  a  third 
light  path. 


The  proper  use  of  this  equipment  required  much  cooperation 
on  the  part  of  the  subject.  Trained  laboratory  personnel  were  there¬ 
fore  used  in  these  experiments.  They  were  also  able  to  apply  their 
knowledge  of  the  method  towards  observations  made  during  the  course  of 
the  experiment  and  thus  give  valuable  information  on  adjustment  of  the 
equipment. 
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Experimental  procedure; 

1)  General  precautions: 

There  is  much  variability  in  the  response  of  the  pupil  from 
one  moment  to  the  next  and  from  run  to  run.  Certain  precautions, 
though,  were  used  to  reduce  the  scatter.  The  experimental  sitting 
was  planned  for  the  morning  after  a  good  night *s  sleep.  The  subject 
was  dark  adapted  U  to  5  minutes  to  a  level  of  illumination  less  than 
the  minimum  which  would  be  presented.  During  the  experiment  room 
light  was  well  controlled  so  that  only  the  stimulating  light  deter¬ 
mined  operating  conditions.  After  the  subject  had  taken  position  for 
an  experimental  run,  he  was  adapted  to  the  conditions  of  the  run  for 
a  period  of  at  least  1J?  secs.  In  the  case  of  the  impulse  experiments, 
two  impulse  trials  were  given  at  the  beginning  of  the  first  trial, 
and  the  run  followed.  Between  runs  the  subject  rested  and  dark 
adapted.  The  atmosphere  was  kept  quiet  to  reduce  psychosensory  phe¬ 
nomenon. 

2 )  Impul se  experiment  s : 

A  function  generator  produced  rectangular  impulses  of 
variable  height  and  width.  A  neutral  density  filter  with  a  transmit¬ 
tance  of  l/lO  was  used  in  the  stimulating  path  and  the  maximum 
stimulating  flux  was  given  an  arbitrary  value  of  one.  Flux  varied 
between  0.1  and  1.0  on  an  arbitrary  scale,  and  impulse  widths  varied 
from  0.01  to  0.U  secs.  Longer  impulse  durations  were  not  used  in 
order  that  a  trial  would  have  as  little  influence  as  possible  on  the 
light  adaptation  of  succeeding  trials. 
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A  trial  was  a  single  application  of  light  at  some  para¬ 
metric  value  and  recurred  in  periods  of  10  secs.  If  there  were  to 
be  five  parameter  changes,  a  five  by  five  latin  square  was  formed  of 
these  values.  A  row  of  this  square  gave  the  succession  of  values 
to  be  used  in  each  run.  There  were  three  replications  at  each  value 
forming  a  group,  and  changes  of  values  were  made  during  the  light- 
off  free  intervals.  In  summary,  then,  for  a  five  by  five  square, 
three  trials  formed  a  group  at  one  value 5  five  groups  at  different 
values  made  a  run,  and  the  experiment  was  composed  of  five  runs. 

The  pen  recorder  was  run  at  a  speed  of  25  ram/sec.  and  gave 
good  resolution  of  the  rising  phase  of  the  response.  Since  the  re¬ 
sponse  of  the  Sandborn  recorder  could  not  reproduce  the  input  pulse 
forms  except  to  indicate  the  foot  of  the  rise  and  the  crown  of  the 
fall,  stimulating  light  output  was  monitored  on  an  oscilloscope. 

In  the  analysis  of  the  data  values  were  read  directly  off 
the  recording  chart  as  relative  pupil  area  (area  x  constant).  Values 
of  maximum  rate  of  rise  were  secured  by  drawing  slopes  to  the  rising 
phase  at  the  point  of  maximum  rate  of  change. 

3)  Sinusoidal  experiments: 

All  experiments  were  run  at  a  frequency  of  l.i*  cps.  This 
was  chosen  as  an  optimum  between  the  irregularities  and  distortions 
seen  at  the  lower  frequencies,  and  the  low  amplitudes  of  response 

found  at  the  higher  frequencies.  Mean  flux  values  were  varied  between 

_2 

extremes  of  transmittance  of  1  to  10  ;  percentages  of  modulation 

ranged  between  0.1  and  0.75.  The  experiment  was  designed  around  a 
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greco-latin  square  in  which  5  values  of  each  parameter  were  randomised 
among  25  runs*  Each  run  was  preceded  by  dark  adaptation  and  was  fol¬ 
lowed  by  a  calibration. 

Analysis  of  the  data  necessitated  subjective  selection  of 
good  pieces.  Particularly  at  lower  frequencies,  a  gaiting  effect  sup¬ 
presses  response,  then  allows  a  period  of  response  to  come  through. 

The  criterion  for  selection  was  to  favor  the  larger  responses  of  the 
run,  but  then  only  provided  that  they  came  in  blocks  of  2-3  waves  and 
had  a  phase  lag  characteristic  for  that  frequency.  The  relative 
amplitudes  were  measured  on  the  recording  graph  and  converted  to  ab¬ 
solute  area  units. 

Statistical  methods: 

The  inherent  biological  variability  of  this  system  required 
the  use  of  replication  and  statistical  methods  to  extract  the  greatest 
amount  of  information  from  the  experiments.  Also,  change  in  the  inter¬ 
nal  operating  conditions  of  the  system  with  time  had  to  be  relegated 
to  an  inferior  position  in  the  determination  of  the  dependencies  under 
investigation.  Towards  these  ends  use  was  made  of  the  experimental 
design  technique  developed  by  R.  A,  Fisher  (1951).  These  techniques 
are  based  on  the  requirement  that  each  parametric  value  be  accorded 
an  equivalent  or  combined  standing  within  the  plan  of  experimental 
grouping  and  coordination.  On  paper  a  design  will  therefore  take  the 
form  of  a  square  array  of  rows  and  columns  of  the  values  of  the  vari¬ 
ables.  Each  value  appears  once  in  each  column  and  each  row  and  is 
assigned  its  place  randomly. 
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In  the  case  of  the  sinusoidal  experiments  a  value  was  two 
dimensional  being  composed  of  percentage  modulation  and  mean  light 
level  values.  The  resulting  array  is  called  a  greco-latin  square. 

To  extract  the  functional  relationship,  rows  and  columns  of  the  array 
were  arranged  to  give  the  values  for  each  parametric  value,  percentage 
modulation  and  mean  level  respectively.  These  were  plotted  to  give 
graphical  representation  to  the  functional  relationships.  Such  data 
can  also  be  analyzed  for  variance  and  thus  a  numerical  comparison  of 
dependence  on  the  two  variables  was  made.  In  the  case  of  the  impulse 
experiments  each  value  of  the  array  was  one  dimensional  and  thus 
there  was  a  replication  of  each  value  by  the  order  of  the  array  (Ijx  or  £x). 
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CHAPTER  VI 


RESULTS 


Impulse  response? 

General  form:  (See  Fig,  5) 

The  rising  phase  is  characterized  by  a  latency,  an  ac¬ 
celerating  phase,  inflection  point,  decelerating  phase,  and  peak. 
Following  the  foot  of  the  stimulus  impulse,  a  0,25  sec,  latency 
intervenes  before  the  initial  rise  of  the  response.  The  function  has 
a  rising  phase  which  reaches  a  maximum  acceleration  in  about  0,2  secs, 
and  then  decelerates  to  the  peak  of  the  response.  Times  to  peak 
averaged  0,56  rfc  .02  secs.,  except  for  responses  to  the  0.1;  sec.  im¬ 
pulse  which  showed  0.1  sec.  longer  duration.  The  falling  phase, 
(redilation)  again  had  a  significant  foot  and  inflection  point. 
Throughout  its  course  it  was  much  slower  than  the  rising  phase.  Later 
segments  of  the  curve  showed  more  variability  in  the  morphology  of 
the  response;  random  behavior  became  more  prominent  and  the  shapes  of 
the  curves  take  on  different  forms.  In  general  the  second  phase 
gradually  decelerated  in  an  exponential -like  manner  shown  in  Fig.  5. 
There  was  a  characteristic  slow  uniform  dilation  in  the  later  phases 
which  makes  curve  fitting  of  the  whole  curve  by  simple  exponential 
functions  unsuitable.  Responses  to  extremes  of  the  input  variables 
did  not  show  large  differences.  Fig.  6. 

Quantitative  data: 

Figure  7  demonstrates  the  variables  in  which  we  will  be 


interested  and  their  notation 
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Fig,  5.  Representative  response  to  pulse.  Pulse  duration 
equals  0.1  sec.;  amplitude  equals  0.5  units.  Abscissa  is 
time  axis  (each  large  division  equals  1/25  sec.);  ordinate  is 
axis  of  amplitudes  (relative). 
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Fig.  6.  Responses  to  extremes  of  pulse  parameters. 
First  two  samples  illustrate  extremes  of  amplitude 
(duration  constant  at  0.05  secs.):  impulse  amplitudes 
equal  0.1  in  (a)  and  1.0  in  (b).  Last  two  samples 
illustrate  extremes  of  duration  (amplitude  constant 
at  0.5  units):  impulse  durations  equal  0.1  in  (c)  and 
0.4  in  (d). 
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Fig.  7.  Impulse  parameters  with  symbolism  used  in  text. 
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Response  height  bears  a  different  relationship  to  impulse 
duration  than  to  impulse  height.  Impulses  of  duration  0.05  or  greater 
show  gradual  increase  of  area  change  with  increasing  duration.  The 
curve  relating  response  to  impulse  height  characteristically  has  a 
small  curvature*  rising  rapidly  from  low  values.  Fig.  8  plots  results 
from  an  experiment  varying  impulse  duration  and  below*  in  Fig.  9  is 
shown  results  from  an  experiment  varying  impulse  height.  A  log  plot 
of  this  data  is  shown  in  Figs.  10  and  11  respectively. 

Maximum  rate  of  rise:  Changing  impulse  width  had  little 
effect  on  the  maximum  rate  of  rise  except  at  the  0.01  sec.  duration 
value,  which  was  significantly  diminished.  There  was,  though*  a  de¬ 
pendence  of  maximum  rate  of  change  on  change  in  impulse  height.  Fig.  12 
plots  rates  of  change  against  response  height  and  the  relationship 
proves  to  be  relatively  linear. 

Results  section  -  Sinusoidal  Experiment" 

A  selected  piece  of  data  from  a  sinusoidal  experiment  is 
illustrated  in  Fig.  13,  and  response  dependence  on  variations  of  input 
percentage  modulation  and  mean  light  levels  are  plotted  in  Figs,  lit  and  . 
15.  The  former  diagram  demonstrates  a  heavy  dependence  of  output  and 
change  on  percentage  modulation.  Dependence  on  mean  light  level  is 
less  significant.  This  observation  is  borne  out  by  variance  analysis 
which  attributes  a  variance  of  1-3  to  variation  in  percentage  modula¬ 
tion,  while  only  0.3  is  due  to  change  in  mean  light  level.  Over 
moderate  ranges  of  operation  the  output  of  the  pupil  system  may  there¬ 
fore  be  said  to  be  relatively  independent  of  mean  operating- level. 
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Fig.  8:  Dependency  on  impulse  duration.  Data  from 
experiment  in  which  impulse  duration  varied;  impulse 
amplitude  constant  at  0.5. 


Fig.  9.  Dependency  on  impulse  amplitude.  Data  from 
experiment  in  which  impulse  amplitude  varied*;  impulse 
duration  constant  at  0.05  secs. 
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Fig.  10.  Log  plot  of  duration  data,  of  Fig.  8. 


Fig.  11.  Log  plot  of  amplitude  data  of  Fig.  9. 
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Fig.  12.  Maximum  rate  of  area  change  vs.  maximum 
response  amplitude.  Relationship  is  relatively  linear. 
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Fig.  13.  Sample  of  sinusoidal  response  data. 
Upper  curve  gives  input  variation  of  light  level 
around  a  mean.  Lower  curve  shows  response. 
Frequency  is  1.4  cps. 
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Fig.  14.  Sinusoidal  amplitude  dependency  on  percentage 
modulation  of  light.  Curve  extrapolated  to  origin. 


Fig.  15.  Sinusoidal  amplitude  dependency  on  mean  light 
level.  Shows  relatively  little  dependence. 


CHAPTER  VII 


DISCUSS  ION 


Sinusoidal  Data: 

The  problem  which  the  sinusoidal  experiments  were  designed 
to  attack  is  the  determination  of  the  nature  of  the  dependency  of 
output  on  the  driving  function.  If  the  output  were  primarily  dependent 
on  the  absolute  amplitude  of  the  input,  one  would  expect  large 
changes  in  output  for  a  ten  fold  change  in  mean  light  level  keeping 
percentage  modulation  constant.  On  the  contrary,  it  has  been  shown 
that  the  output  did  not  increase  with  an  increase  in  mean  level,  but 
in  fact  probably  fell  off  at  higher  levels.  On  the  other  hand  a 
change  in  the  percentage  modulation  produced  a  corresponding  change  in 
the  same  direction  in  the  output.  Since  zero  percent  modulation  would 
produce  zero  output,  the  output  vs.  input  percent  modulation  curve 
shown  in  Fig.  lit,  p.  1*0,  crosses  the  positive  absissa.  For  purposes 
of  curve  fitting  we  have  taken  the  logical  step  of  drawing  the  curve 
through  the  origin.  It  will  be  noted  that  this  relationship  is  non¬ 
linear. 

In  summary,  then,  the  significance  of  the  results  of  these 
experiments  may  be  simply  stated:  Within  the  range  investigated,  the 
significant  driving  function  is  not  the  absolute  amplitude  but  rather 
the  percentage  amplitude  change  of  the  input  sinusoid. 

Pulse  data:  general  discussion 

Since  the  pupil  system  has  already  been  considered  from  the 
standpoint  of  linear  analysis  with  regard  to  response  in  the  frequency 
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domain,  it  would  be  well  to  begin  our  evaluation  of  the  pulse  ex¬ 
periments  in  terms  of  comparable  linear  systems  suggested  by  those 
investigations.  The  weight  of  this  suggestion  is  pressed  upon  us  be¬ 
cause  linear  theory  demands  a  one  to  one  correspondence  between  the 
two  domains,  i.e,  suitable  time  data  and  frequency  data  will  both 
completely  describe  a  system  and  results  from  one  method  transform 
into  that  of  the  other.  Stark  and  Sherman  suggested  that  their  data 
might  be  interpreted  as  indicating  a  three  lag  system,  and  for  their 
analysis  they  called  each  lag  approximately  equal  to  0.1  sec.  We 
will  therefore  discuss  the  response  of  two,  three,  and  four  lag 
linear  systems  to  single  pulse  inputs. 

Systems  which  may  be  represented  by  first  order  differential 
equations,  or  single  lag  systems,  respond  to  a  pulse  by  an  exponential 
rise  through  the  duration  of  the  pulse  with  an  exponential  fall  im¬ 
mediately  following.  The  data  on  the  present  system  had  a  peak  of 
response  well  after  the  cessation  of  the  impulse,  and  there  is  also  a 
significant  foot  in  the  rising  phase.  Single  lag  systems  can  therefore 
be  excluded  from  consideration. 

Curves  in  time  of  multilag  systems  show  a  gross  similarity 
but  differs  on  closer  inspection.  Increasing  order  increases  the 
extent  of  the  foot  of  the  rising  phase  and  the  rapidity  of  descent  of 
the  falling  phase  for  a  given  duration  to  peak.  These  characteristic 
differences  are  illustrated  in  the  superimposed  plots  of  two,  three 
and  four  lag  systems  with  equal  time  constants  shown  in  Fig.  16.  Time 
constants  and  gains  have  been  adjusted  to  give  equal  peaks  of  response. 
The  scale  of  the  coordinates  is  arbitrary. 
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Fig.  16.  Pulse  responses  of  ideal  2,  3,  and  4  lag 
systems.  Gain  factor  and  time  constant  adjusted  to 
make  peaks  coincide. 


Fig.  17.  Representative  pupil  response  compared  with 
ideal  lag  systems.  Input  pulse  duration  was  0.05  secs.; 
amplitude  was  0.2  units. 


-  hh  - 


In  Fig,  17  a  reproduction  of  a  response  of  the  pupil 
system  is  compared  with  the  ideal  responses  of  Fig,  16.  Peak 
responses  have  been  adjusted  to  coincide.  It  is  evident  that  the 
pupil  system  closely  resembles  a  third  order  system,  in  excellent 
correspondence  with  the  hypothesis  derived  from,  frequency  domain 
studies. 

Critical  quantitative  characteristics: 

A  fundamental  notion  in  linear  analysis  is  the  impulse  re¬ 
sponse  of  a  system.  The  impulse  input  is  a  pulse,  the  duration  of 
which  has  been  reduced  to  an  infinitesimal  interval  and  the  magnitude 
increased  infinitely  in  such  a  way  that  the  energy  content  of  the 
function  equals  one.  In  practice  the  impulse  function  may  be  ap¬ 
proximated  by  pulses  with  durations  less  by  an  order  of  magnitude 
than  the  time  constant  of  the  system.  Magnitude  of  response  of  linear 
systems  to  such  inputs  is  proportional  to  energy  content  (amplitude  x 
duration)  of  the  pulse  and  response  duration  to  peak  is  a  minimum  for 
a  given  system.  Pulses  of  longer  duration  produce  responses  which  fall 
away  from  the  energy-magnitude  proportionality  and  in  which  the  duration 
increases  beyond  that  seen  for  impulses.  For  long  pulse  durations  the 
response  amplitude  asymptotically  approaches  that  seen  in  responses  to 
a  step  while  the  duration  to  peak  approaches  the  duration  of  the  pulse. 
These  relationships  are  diagrammed  schematically  for  a  multilag  system 
in  Figs.  18  and  19. 

Impulse  time  to  peak  is  closely  related  to  the  time  constants 
of  the  system.  In  a  system  with  equal  lags  the  relationship  is  simply 
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Fig.  18.  Maximum  amplitude  vs .  input  pulse  duration. 
Diagramatic  representation  of  arbitrary  physical  system. 
Curve  asymptotically  approaches  step  response  amplitude. 


Fig.  19.  Duration  to  peak  response  vs.  input  pulse  duration. 
Diagramatic  representation  of  arbitrary  physical  system. 
Curve  asymptotically  approaches  duration  of  pulse. 
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the  proportionality 


where  n  is  the  order  of  the  system. 

Comparison  of  the  pupil  system  to  two,  three  and  four  lag 
systems  is  illustrated  in  Figs,  20  and  21,  in  which  experimental  points 
are  represented  by  circles.  Maximum  response  height  increases  only 
gradually  on  a  plateau  with  increasing  duration  and  only  at  the  shortest 
duration  pulse  is  there  an  indication  of  a  break  from  the  high  main¬ 
tained  value.  Experimental  values  for  duration  to  peak  fall  below 
the  curve  which  represents  this  relationship  in  the  linear  systems. 

It  is  important  to  note  that  the  extrapolated  value  of 

f 

the  durations  to  peak  is  about  0.50  secs.  If  the  pupil  system  had 
three  equal  lags,  this  would  indicate  that  the  values  of  the  time 
constants  would  be  0.2$  secs. 

Response  as  a  function  of  input  pulse  amplitude 

Fig.  9 ,  p.  36,  indicates  that  amplitude  is  a  non-linear 
function  of  the  input.  In  accordance  with  the  Weber-Fechner  law 
response  of  sensory  systems  frequently  is  compared  to  the  log  of  the 
input  amplitude.  As  seen  in  the  log  plot  of  the  data  on  this  system 
illustrated  in  Fig,  11,  the  relationship  is  still  not  a  proportionality. 
As  was  mentioned  in  Chapter  V  [Methods],  though,  there  was  a  small 
minimum  light  level  present  on  which  the  pulses  were  superimposed,  and 
this  might  have  changed  the  level  of  retinal  adaptation.  Response  to 
small  pulses  would  be  affected  more  than  that  to  large  pulses  resulting 
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Fig.  20.  Amplitude  vs.  pulse  duration  in  pupil  and  ideal 
2,  3,  and  4  lag  systems.  Gain  of  pupil  system  arbitrary. 
Time  constants  of  lag  systems  are  equal  and  values  obtained 
from  duration  to  peak  data,  Fig.  21. 


Fig.  21.  Duration  to  peak  vs.  pulse  duration  in  pupil  and 
ideal  2,  3,  and  4  lag  systems.  Pupil  system  represented  by 
lower  curve  through  data.  Time  constants  of  ideal  systems  are 
equal  and  adjusted  to  give  impulse  duration  equal  to  that  of 
the  pupil  system. 


in  lower  amplitudes  at  those  values.  Thus,  system  response  might 
well  be  proportional  to  the  log  input  in  this  range  of  operation. 

An  argument  against  an  hypothesis  of  mechanical  range 
limiting  in  the  motor  system  producing  amplitude  saturation,  similar 
to  that  used  by  van  der  Tweel,  is  based  on  the  maximum  rate  of  response 
vs.  amplitude  of  response  relationship.  If  there  were  amplitude  satur¬ 
ation,  one  would  expect  a  maximum  response  rate  vs.  amplitude  curve 
as  illustrated  in  Fig.  22.  This  is  to  be  compared  with  the  relation¬ 
ship  in  a  velocity  saturation  system  illustrated  in  Fig.  23.  The 
curve  of  experimental  values  shown  in  Fig.  12,  p.  38,  clearly  does 
not  have  an  upward  concavity,  and  in  fact  might  have  a  slight  downward 
concavity  suggesting  velocity  saturation.  This  is  of  interest,  since 
Stark  (1961)  has  postulated  velocity  saturation  in  this  system  to 
account  for  data  derived  in  a  totally  different  way.  The  evidence  is 
therefore  against  mechanical  range  limiting  in  this  operating  range. 

Critique  of  the  linear  lag  hypothesis 

The  form  of  the  response  of  the  pupil  system  to  pulses  fits 
very  well  with  the  hypothesis  that  the  system  consists  of  three  equal 
lag  elements  as  the  primary  response  shaping  parameters.  A  log-like 
amplitude  operator  could  well  be  incorporated  into  such  a  scheme  as 
long  as  it  came  earlier  in  the  system  than  the  frequency  dependent  ele¬ 
ments,  required  by  the  argument  presented  by  Selin.  The  present 
evidence  would  indicate  that  the  time  constants  were  equal  to  0.25  secs. 
This  value  lies  at  the  upper  extreme  to  the  range  for  the  time  constant 
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Fig.  22.  Maximum  rate  vs.  maximum  amplitude 
relationship  in  system  with  mechanical  limit 
saturation. 


Fig.  23.  Maximum  rate  vs.  maximum  amplitude 
relationship  in  system  with  velocity  saturation. 
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suggested  by  Stark  and  Sherman,  but  it  would  be  an  even  more  suitable 
value  for  our  own  amplitude  vs.  frequency  data.  From  the  work  of 
Stark  and  Baker  (1958),  who  demonstrated  a  large  shift  in  frequency 
of  natural  oscillation  following  the  application  of  drugs  to  the 
cornea,  the  proposed  system  would  most  likely  have  the  major  lag 
components  at  the  periphery. 

The  data  on  the  maximum,  amplitude  and  the  duration  to  peak 
as  functions  of  pulse  duration,  though,  would  indicate  that  other 
processes  are  in  operation  not  accounted  for  by  this  hypothesis.  While 
it  is  not  necessary  to  discard  the  above  concepts  in  the  light  of  this 
evidence,  certain  modifications  are  required  to  make  the  theory  more 
complete.  Although  the  present  experimental  work  presents  insufficient 
data  to  justify  construction  of  an  elaborate  theory,  the  remainder  of 
this  chapter  will  discuss  what  might  be  a  fruitful  point  of  view 
towards  this  system,  accounting  for  the  latter  data  and  yet  preserve 
the  afore  mentioned  hypothesis. 

A  proposed  mechanism: 

If  energy  is  taken  to  be  the  significant  driving  function, 
the  experimentally  derived  relationship  of  amplitude  vs.  pulse  duration 
makes  evident  that  short  pulse  durations  are  favored  over  long  durations 
in  the  response  of  the  system.  Such  relationships  have  been  well 
studied  and  documented  in  the  physiology  of  sensory  elements  (See  the 
Handbook  of  Physiology,  section  of  neurophysiology,  V.  I).  Classically, 
a  receptor  system  initially  responds  to  a  stimulus  with  high  frequency 
firing  which  decreases  with  time  to  either  a  low  maintained  value  or  to 
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zero.  Hartline  and  Grabm  (1932)  showed  this  process  to  take  place 
in  single  receptor  elements  in  the  arthropod  eye.  Riggs  et  al  (1953) 
have  demonstrated  the  adaptation  process  in  the  human  eye  by  main¬ 
taining  the  projection  of  a  stimulus  on  one  set  of  receptors  using  an 
optical  arrangement.  Hartline  et  al  (1959*  review)  have  analyzed  a 
process  of  adaptation  in  the  eye  of  Limulus  which  is  accounted  for 
solely  by  neuronal  interraction  mechanisms.  Similar  physiological 
arrangements  exist  in  the  mammalian  eye  which  could  result  in  similar 
properties. 

In  the  light  of  the  experimental  evidence  and  these  physio¬ 
logical  mechanism s 3  which  might  well  be  operating  in  this  system,  one 
can  conceive  of  an  operator  early  in  the  system  such  as  that  depicted 
in  Fig.  2lu  Short  duration  impulses  would  be  to  a  large  part  passed 
but  the  later  phases  of  long  duration  pulses  would  be  greatly  at¬ 
tenuated.  By  proper  adjustment  of  the  parameters,  the  convolution  of 
such  a  function  with  the  proposed  lag  system  could  duplicate  the 
characteristics  of  maximum  amplitude  vs.  pulse  duration  curve.  The 
presence  of  such  an  operator  would  also  explain  the  experimental 
duration  to  peak  vs.  pulse  duration  curve.  With  such  an  operator,  in 
principle  well  within  the  confines  of  neurophysiological  mechanisms, 
combined  with  the  concept  of  a  system  of  three  equal  lags,  the  data 
here  presented  can  be  well  explained. 
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Fig.  24.  Possible  input  operator.  General  form 
of  hypothesized  operator  which  could  account  for 
pulse  data. 


CHAPTER  VIII 


CONCLUSIONS 

The  problem  in  question  has  been  to  investigate  the  func¬ 
tional  dependencies  of  pupil  area  on  light  flux  at  the  plane  of  the 
retina.  Experiments  in  two  different  modes  of  system  operation  have 
given  us  separate  pieces  of  information  concerning  this  relationship. 
Results  of  the  sinusoidal  experiments  have  not  been  theoretically  ex¬ 
plored.  In  the  discussion  of  the  pulse  experiments  the  view  point 
that  the  response  of  the  system  is  primarily  shaped  by  a  multilag 
network  is  modified  to  generate  a  relationship  more  consistent  with 
the  present  data.  While  the  idea  of  the  time  dependent  input  function 
accords  with  neurophysiological  concepts,  no  attempt  has  been  made  to 
formalize  this  relationship.  The  present  study  only  hints  at  the 
form  of  the  function  and  does  not  offer  enough  data  for  a  good  analysis. 
For  these  reasons  the  hypothesis  has  only  been  touched  upon,  and  no 
further  investigation  has  been  made  of  the  consequences. 

Specifically,  an  attempt  has  not  been  made  to  correlate  the 
pulse  response  with  amplitude  dependent  sinusoidal  response  character¬ 
istics.  Yet,  as  previously  mentioned,  Tyler  has  pointed  out  that  the 
rate  of  rise  of  the  response  to  a  step  function  is  not  consistent  with 
the  limitations  implied  by  the  frequency  response  characteristics  of 
this  system.  It  is  probable,  then,  that  the  system  is  operating  in 
different  states  under  these  two  different  modes  of  stimulation.  Also, 
Stark  has  demonstrated  that  there  is  hardly  a  break  in  the  amplitude  of 
response  to  a  sine  wave  input  when  the  mean  level  is  cut  by  an  order  of 
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magnitude.  This  implies  that  the  accommodation  process  accounting 
for  this  phenomenon  is  rapid.  These  two  observations  make  reasonable 
an  identification  of  the  dependency  of  sinusoidal  output  on  percentage 
light  change  and  the  hypothe sized  function  in  the  impulse  analysis. 

Much  more  work,  though,  needs  to  be  done  with  the  pulse  experiments 
before  such  a  relationship  can  be  established.  In  particular,  a 
greater  range  of  pulse  durations  must  be  explored,  as  well  as  the 
effect  of  light  adaptation  on  pulse  response. 

In  the  final  analysis  the  goal  of  these  studies  is  to  obtain 
the  complete  functional  relationship  between  light  in  and  area  out. 
Input-output  analysis,  along  with  ideas  of  basic  physiological  mechanisms, 
has  lead  to  an  increased  understanding  of  these  relationships,  or 
at  least  has  presented  some  rational  formulations.  It  is  not  to  be 
expected  that  systems  analysis  alone  will  lay  open  the  fundamental 
mechanisms  of  the  pupil  system.  Indeed,  such  techniques  are  geared 
rather  to  illucidate  functional  forms  of  response  and  in  this  capacity 
tests  hypotheses  and  makes  suggestion  as  to  the  nature  of  the  funda¬ 
mental  mechanisms.  Thus  used,  this  method  has  much  more  to  offer  in 
assisting  the  investigation  of  complicated  signal  flow  systems. 
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CHAPTER  IX 


SUMMARY 

The  advantages  of  using  the  pupil  system  as  a  model  for  a  study 
of  neurological  reflex  are  presented. 

Relevant  anatomy  and  recent  physiological  investigation  is  re¬ 
viewed. 

The  point  of  view  of  systems  analysis  is  developed*  and  illustration 
made  of  use  to  which  it  has  been  put  in.  study  of  the  pupil  system. 
The  experimental  methods  of  the  present  study  are  reviewed*  and 
designs  for  experiments  clarifying  the  dependency  of  the  response 
on  input  parameters  for  the  sinusoidal  and  pulse  modes  of  stimu¬ 
lation  are  outlined. 

Results  of  the  sinusoidal  experiments  indicate  that  amplitude  of 
response  is  primarily  dependent  on  percentage  light  change  rather 
than  the  absolute  value  of  light  change. 

Results  of  the  pulse  experiments  demonstrate  that  the  form  of  the 
response  fits  very  well  the  time  curve  of  response  of  a  third 
order  physical  system. 

The  dependency  of  response  on  pulse  amplitude  indicates  that  there 
is  probably  a  log  operator  early  in  the  signal  flow  path.  A 
linear  maximum  rate  vs.  maximum  amplitude  relationship  supports 
the  hypothesis  that  there  is  minimal  mechanical  limiting  of  the 
iris  in  this  range  of  operation. 

The  characteristics  of  second*  third  and  fourth  order  linear 
systems  are  reviewed*  and  it  is  shown  that  the  pupil  system 
cannot  be  directly  compared  to  a  network  of  simple  lags.  An 
hypothesis  is  presented  to  account  for  this  discrepancy. 
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